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Abstract. In low-mass Red Giant Branch (RGB) and Asymptotic Giant Branch (AGB) stars, 
anomalous mixing must transport material near the hydrogen-burning shell to the convective en- 
velope. Recently, it was suggested that buoyant magnetic flux tubes could supply the necessary 
transport rate [1]. The fields are assumed to originate from a dynamo operating in the stellar inte- 
rior. Here, we show what is required of an a — £2 dynamo in the envelope of an AGB star to maintain 
these fields. Differential rotation and rotation drain via turbulent dissipation and Poynting flux, so if 
shear can be resupplied by convection, then large-scale toroidal field strengths of (Bq) ~ 3 x 10 4 G 
can be sustained at the base of the convection zone. 
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INTRODUCTION 

In low-mass Red Giant Branch (RGB) and Asymptotic Giant Branch (AGB) stars, an 
anomalous transport mechanism must cycle material from the convective envelope to the 
stable radiative region where nuclear processing can occur |P,|3Q. In particular, material 
must be pumped downward from the convection zone, exposed to partial H burning 
(so-called Cool Bottom Processing; CBP) and returned to the convective envelope. 
Isotopic abundances provide constraints on the required processing temperatures and 
mass transfer rates. Recently, it was proposed that buoyant magnetic flux tubes can 
provide the necessary transport yj]. 

In evolved stars, magnetically mediated outflows have been proposed as the origin of 
bipolaritv in post-Asympotitc Giant Branch stars (post- AGB) and planetary nebula (PN) 
HBfiB]- The needed large-scale field likely results from a dynamo operating during 
the RGB or AGB. In both phases, a convective envelope coupled with differential rota- 
tion in the stellar interior can amplify the field via an a — Q. dynamo QV[] . However, as the 
field amplifies, differential rotation energy is drained and rapidly terminates the dynamo 
(< 100 yrs) yLyQ. This would favor binary shaping mechanisms for post-AGB/PNe and 
a binary induced dynamo (see O, [80 for detail). Only if convection resupplies differen- 
tial rotation (analogous to the A -effect in the Sun [Q]) could a sustained dynamo operate 
throughout the RGB or AGB phases without a binary [7]. While sustained fields are re- 
quired to shape post-AGB outflows, the dynamo required for mixing may be weaker and 
need not influence the outflow. It may be that a weaker dynamo (supplying sufficient 
mixing) is sustained throughout the AGB phase, but that a binary companion is required 



1 nordhaus@pas.rochester.edu 



FIGURE 1. Figure from [7]. A meridional slice of dynamo geometry. The left figure shows the global 
geometry of the AGB star. The right figure is a close-up view of the dashed region on the left. The en- 
effect is driven by convection and occurs in layer of thickness L\ above the differential rotation zone. 
The poloidal component of the field is pumped downwards into the differential zone, where it is wrapped 
torodially due to the £2-effect. 



to power post- AGB bipolarity. 

Here, we present a dynamical a — Q. dynamo operating at the base of the convection 
zone in a 3 M & AGB star. We consider two models: (i.) shear is not resupplied, and a 
transient dynamo results (ii.) convection resupplies shear, sustaining the dynamo. For 
both models, the field penetrates the shear zone and may transport material to the H 
burning shell. We compare our results to the field strengths and penetration depths 
needed for viable magnetic mixing in low-mass RGB and AGB stars [1]. We comment 
on future directions. 

AN EVOLVED STAR DYNAMO 

A 2-D schematic of our 1-D, 3 M & AGB statH is presented in Fig. [TJ In such a star, the 
convection zone extends from the stellar surface to the interface between the convective 
and radiative zones. In this layer, convective twisting motions convert buoyant toroidal 
field into poloidal field through the oc-effect. Below the convection layer, the differential 
rotation zone extends to the core and shears poloidal field back into toroidal field through 
the ^-effect. 

The differential rotation profile in the AGB interior is unknown. In order to employ 
a reasonable initial profile, we assume that the main sequence surface rotation velocity 
is 200 km s . A radial shear profile is then obtained by conserving angular momentum 
on spherical shells during post-MS evolution Q6Q . 

We solve for the time evolution of dynamic quantities at the interface between the 
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convective region and shear zone (r = r c ). The rotation profile across the shear layer 
varies from Q. at the interface to £1 + A£l at r = r c —L. Thus, A£l is a measure of shear 
in the differential rotation zone. If AO. = 0, the system exhibits solid body rotation. 
The average poloidal field, {Bp), and average toroidal field, ufy) amplify from 1 G 
seed values until they are quenched through a drain of the available differential rotation 
energy. We refer to reader to QVU for a derivation of the precise 1-D, mean-field equations 
solved. 



Shear Zone Penetration Depth 

To capture aspects of 2-D, mean-field geometry within the framework of our simple, 
1-D time-dependent model, we employ two turbulent diffusion coefficients: fi p , corre- 
sponding to diffusion of the poloidal field (which grows primarily in the convective 
region) and j8^, corresponding to diffusion of the toroidal field (which is amplified in 
the differential rotation zone; see Fig. 1). We also use as the turbulent diffusion co- 
efficient for the toroidal velocity. Since the convective region is highly turbulent and the 
differential rotation zone is more weakly turbulent, we take jSa <C f5 p . The value of jS* 
governs how far the poloidal component of the field can diffuse into the shear zone in a 
cycle period. The further into the shear zone the poloidal field can penetrate, the greater 
the shear energy that can be extracted and utilized by the dynamo. The coefficient, J[L 
also governs how much heat is generated by turbulent dissipation in the shear layer. 

To quantitfy this, we define 8 as the depth at which the poloidal field can diffuse into 
the shear layer in a cycle period, i. We have 

5 ~ 06^)1 (1) 

If 8/L = 1 , the poloidal component of the field reaches the core boundary and extracts all 
available shear energy. If 8/L < 1, then the poloidal field extracts a fraction of the total 
shear energy and can not penetrate as deeply. Thus for a given /3^, we can determine the 
depth to which the field diffuses. As an aside, we note that for the maximally growing 
dynamo mode, % is weakly dependent on and hence, would change the power index 
slightly. Although we do not expect this effect to be significant, future work should use 
the maximally growing mode. 



The Need For Convection Resupplying Shear 

To investigate the back-reaction of field growth on the shear, we apply our dynamical 
model using the initial rotation profile previously described. Our results are shown in 
Fig. [2l The dynamo terminates rapidly as the shear energy is drained after ~ 20 yrs. 
In addition, the loss of Poynting flux spins down the interface by a factor of ~ 3. In 
order to sustain the dynamo through an AGB lifetime, a constant differential rotation 
profile must be established. This occurs in the sun as convection re-seeds shear through 
the A -effect [9]. Although it remains to be established if a similar effect occurs in 
evolved stars, by analogy to the solar case, we allow a fraction of the turbulent energy 
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FIGURE 2. Figure from The differential rotation energy is allowed to drain through field ampli- 
fication and turbulent dissipation. We define, B^pp^is] = JoE^ PFdis ^ (t')dt' and label M (Epf) and T (Edis) 
on the top right plot to distinguish between the thermal and magnetic contributions. If £[pF,dis]/Ebmd > 1, 
then the AGB envelope was unbound via heat from turbulent dissipation or time-integrated Poynting flux. 
Peak field strengths are a factor of ^ 5 — 10 less then those obtained in (05|]). Differential rotation energy 
is drained in < 20 yrs. Lowering the turbulent dissipation coefficient in the shear region results in the 
differential rotation energy draining at a slower rate, allowing the field to sustain for longer periods of 
time (~ 40 — 50 yrs). However, peak field strengths remain the same and are a function of AQ, and Q.. 



cascade to resupply shear. Additionally, we keep the rotation at the interface fixed. This 
is physically equivalent to storing Poynting flux in the interface region. If the field is 
trapped, Poynting flux does not emerge from the layer and thus, does not spin down the 
envelope. An AGB dynamo will thus sustain when the following two conditions are met: 
(i.) convection resupplies shear (ii.) Poynting flux is stored inside the envelope tn\- 

For our 3 M Q AGB star, under those conditions, a steady-state dynamo is established 
when ~ 1% of the turbulent cascade energy per unit time reinforces the shear (see Fig. 

REQUIREMENTS FOR MAGNETIC MIXING 

For a given bulk transport rate, Cool Bottom Processing predicts the chemical and 
isotopic evolution of the convective envelope. Without appealing to a specific driving 
mechanism, CBP yields explicit temperature requirements and nuclear processing times 
needed to match observations for a specific stellar model. To investigate the physical 
origin of CBP, hydrodynamical models, such as rotationally induced meridional mixing 
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FIGURE 3. Figure from 10]. Convective resupply results in a steady-state differential rotation profile. 
For the left column, the envelope of the poloidal, toroidal and Poynting flux is plotted. The Poynting 
flux is sustained at ~ 5 x 10 34 erg/s. The sustained Poynting flux supplies enough energy to unbind the 
envelope of our 3 M© model at the end of the AGB phase (~ 10 5 yrs). In this figure, ~ 0. 1 % of the cascade 
power must be used to sustain the differential rotation to supply the requisite Poynting flux. This choice 
of parameters corresponds to a magnetically dominated explosion at the end of the AGB. 



were suggested. However, these models do not supply sufficient mixing to explain 
observations of low-mass RGB field and globular cluster stars ifioll . 

Recently, as an alternative to hydrodynamic models, magnetically induced mixing 
was proposed [1]. Buoyant toroidal fields in the radiative zone carry material to the 
convective envelope and drive mass circulation. By requiring CBP transport rates, it 
is possible to constrain magnetic field strengths and buoyant velocities necessary for a 
given stellar model. 

For an initial 1.5 M & star with metallicity one-half solar [[Till , constraints are cal- 
culated for the subsequent RGB and AGB phases. In both phases, subadiabatic zones 
similar to the solar tachocline exist J12Q. Since our previous dynamo calculations apply 
to an AGB star, we focus on that phase. 

For the AGB phase, the total mass has decreased to M = 1.2 M Q . The base of the 
convection zone is located at r c = 5.7 x 10 10 cm. For a mass flux of M = 10~ 6 M & /yr, 
material must mix to a radius of r C BP = 2 x 10 9 cm which is slightly above the location 
of the H burning shell (r H = 1.5 x 10 9 cm) |Q]]. 

Given a mass transfer rate, mixing depth and buoyant rise velocity, the magnetic field 
strength at r# and r c can be constrained. If the field is purely toroidal and in flux tubes 
that rise individually, then B t — B (rcBp) = 5 x 10 G yj]. At the base of the convection 



zone, this corresponds to a field of B c = B (r c ) ~ 9 x 10 3 G for a dipole geometry. 

Direct comparison to the dynamo calculation is difficult as our AGB star is 2.5 times 
more massive. In addition, our 3 M & structure corresponds to a star that has just reached 
the AGB. Future work should investigate a dynamo operating in a lower mass star, during 
the interpulse periods. Nevertheless, for plausibility arguments, a general comparison 
with our sustained dynamo is useful. 

From our 3 M Q AGB star, the toroidal field strength at the base of the convection 
zone is (B^ jC ) ~ 3 x 10 4 G, suggesting that it may be possible to produce sufficiently 
high average fields. However, since we do not solve for the radial structure, we cannot 
determine the field strength where CBP nucleosynthesis occurs. 

We can, however, estimate the depth to which the magnetic field penetrates the shear 
layer. Based on the initial conditions used in Fig. [3] Eqn. 1 gives that the field reaches a 
radius of ~ 3 x 10 10 cm with approximately 15 cycles per year. This is encouraging and 
suggests that significant penetration into the shear zone may occur. 



Flux Tubes versus Average Fields 

Our dynamo calculations are mean-field and correspond to volume averages. When 

(p) 

the ratio of average thermal to magnetic pressure /3 = jjpxrj^ >> 1 it is possible for the 

magnetic field to be concentrated in flux tubes with an overall volume filling fraction of 

l 



the magnetic field f « 1, to be small, such that / ~ i 111311 . When flux tubes are in 



presure balance with their exterior, the mean-field strength is then related to a flux tube 



field, B t by (B 2 ^ ~ / (Bf} fl 1 3M - Since the field constraints on magnetic mixing from 
Ref. 01] constrain the magnitude of flux tube fields with total / << 1, the mean field 
can be much lower than the flux tube field, relaxing the demands on the strength of the 
dynamo. 

Further investigation of all of the above is warranted for relevant stellar models. 
Determining the radial diffusion coefficients will constrain the amount of shear energy 
available to the dynamo and thus, determine the penetration depth and field strengths. 
For a given stellar model, constraints on the fraction of convective energy needed to 
sustain a steady shear profile should be calculated in both the RGB and AGB phases. In 
the RGB phase, weaker fields are required. Realistic rotation profiles should be utilized 
when known. Finally, > 2-D model calculations incorporating Q-quenching will provide 
radial information on magnetic field strengths. 



CONCLUSIONS 

It was recently suggested that buoyant toroidal fields, amplified in the stellar interior, 
may supply the necessary chemical transport rates from the hydrogen burning shell to 
the convective envelope required from Cool Bottom Processing in low-mass RGB and 
AGB stars H. 

In this paper, we have presented results of a dynamic a — Q. dynamo operating in 
an AGB star. Amplification of the large-scale field drains shear energy rapidly (< 100 



yrs.). For an isolated star, stringent conditions (re-supply of shear + storage of Poynting 
flux) must be met if a dynamo is to sustain through the RGB or AGB phases. This can 
be accomplished if ~ 1% of the turbulent energy cascade reinforces shear. In this case, 
our AGB model produces strong toroidal fields (B<p ~ 3 x 10 4 G) at the base of the 
convection zone that can diffuse into the shear layer. This suggests that magnetic fields 
may be a plausible mechanism for mixing in low-mass RGB and AGB stars, however 
future work is needed before this can be firmly established. 
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